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We present a shortened and simplified version of our proof [f| of the uniqueness of the scaling 
solution for the infrared asymptotics of Green functions in Landau gauge Yang-Mills theory. The 
simplification relates to a new RG-invariant arrangement of Green functions applicable to general 
theories. As before the proof relies on the necessary consistency between Dyson-Schwinger equations 
(DSEs) and functional renormalisation group equations (FRGs). We also demonstrate the existence 
of a specific scaling solution for both, DSEs and FRGs, that displays uniform and soft kinematic 
singularities. 
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I. INTRODUCTION 

The infrared behaviour of strongly interacting quan- 
tum field theories is of general interest. In particular the 
infrared behaviour of Landau gauge Yang-Mills theory is 
a much debated issue in the past years. This is due to the 
close connection of the ghost and gluon propagators to 
confinement scenarios proposed by Kugo and Ojima [l|, 
Gribov Q and Zwanziger 0]. Within functional meth- 
ods a self-consistent infrared asymptotic solution of the 
whole tower of Dyson-Schwinger equations (DSEs) and 
functional renormalisation group equations (FRGs) has 
been found [3, 0, @, 0, B, Q that supports these scenar- 
ios and is also consistent with global BRST symmetry 
[lfj. This scaling solution implies that all Greens func- 
tions scale like a power of momentum in the infrared 
with interrelated anomalous dimensions if all momenta 
are scaled. That entails the absence of mass scales be- 
low which some degrees of freedom decouple. In such 
a scenario the scaling power of Green functions can be 
extracted by a power counting analysis. As a further di- 
rect consequence all couplings have fixed points at zero 
momentum. 

An alternative infrared solution of Yang-Mills th eory 
is provided by the decoupling solution [ll|, [l2|, HH, (Hf. 
This type of solutions has been discussed in detail in 
[Io| . and has been shown to be inconsistent with global 
BRST symmetry. In the present context it can be in- 
terpreted as the decoupling of (part of) the propagating 
degrees of freedom below a mass scale. In such a case 
the infrared asymptotics cannot be fixed uniquely by a 
scaling analysis. 

In 0] we suggested a combined analysis of the towers 
of DSEs and FRGs for an infrared scaling analysis, be- 
ing applicable to general theories. We have shown that 
apart from decoupling there is only one, unique scaling 
solution of infrared Yang-Mills theory in Landau gauge. 
Here we present a greatly simplified version of our proof 



which also allows the reduction of the number of pre- 
suppositions. The result, of course, is the same as in 
0, B @ • In addition we present an explicit scaling solu- 
tion for DSEs and FRGs that involves also kinematical 
singularities [l^]. The knowledge of these kinematical 
singularities is not necessary for the proof of the exis- 
tence and uniqueness of the global scaling, and were not 
discussed explicitly in our previous work y|. 

In Section|TT]we introduce FRG and DSE equations for 
the effective action and discuss momentum and RG scal- 
ing for the effective action and its vertices. With help of 
Appendix[A]the functional equations are written in a sim- 
ilar form. A convenient parameterisation of the vertices 
is introduced that splits off the renormalisation group 
(RG)-scaling and reduces the discussion to that of the 
scaling properties of RG-invariant quantities. This very 
natural reduction is the key ingredient of the simplifica- 
tion of the proof, and is applicable to general theories. It 
also allows us to provide heuristic arguments for the ex- 
istence of a unique scaling solution, that should facilitate 
the following of the proof. In Section IIIII we derive the 
respective scaling constraints from FRG and DSE equa- 
tions, the combination of which provides a unique scaling 
solution. In Section IIVI we extend our analysis to kine- 
matical singularities of vertices, the details can be found 
in Appendix [B] We close with a short summary of our 
findings. 



II. FUNCTIONAL RELATIONS FOR THE 
EFFECTIVE ACTION 

The starting point of our analysis is the functional form 
of FRGs and DSEs, depicted in Fig. Q] and Fig. [U re- 
spectively. In Fig. [1] we have rewritten the FRG in 
a form similar to the DSEs in Fig. [2l This is detailed 
in Appendix [X] It is well-known from the evaluation of 
critical physics that the DSE is less amiable to the dis- 
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FIG. 1: Infrared asymptotics of the FRG. Filled circles de- 
note fully dressed field dependent propagators. Empty circles 
denote fully dressed field dependent vertices. 
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FIG. 2: Functional Dyson-Schwinger equation (DSE) for the 
effective action. Filled circles denote fully dressed field de- 
pendent propagators. Empty circles denote fully dressed field 
dependent vertices, dots denote field dependent bare vertices. 



cussion of scaling than the FRG. In general scaling in the 
DSE only comes from a combination of diagrams which 
hosts cancellations effectively leading to the substitution 
of the bare vertex present in each diagram with dressed 
ones. In general such a cancellation also includes the 
classical term. One example for such cancellations is the 
</> 4 -theory. 

In the present case, however, the functional ghost DSE 
(second line of Fig. [2]) prohibits cancellations between di- 
agrams as it only consists of one diagram. It is for this 
reason that the DSE system of Yang-Mills theory can- 
not be subject to cancellations except for kinematical 
and global symmetry reasons. Moreover, in both sets 
of equations Fig. [TJ Fig. [5] we have the classical term in 
the DSE or the initial condition for the FRG 5S(a)/54>- 
For infrared enhanced vertices and inverse propagators 
this term is subleading. On the other hand, scaling of 
infrared suppressed vertices and inverse propagators re- 
quires cancellations on the right hand side of the FRGs 
and DSEs between the diagrams and the classical term. 
As will become clear later the only place where such a 
cancellation necessarily has to occur for a scaling solu- 
tion are the DSE and FRG of the ghost propagator. The 
presence or absence of cancellations in these equations 
therefore decides about the existence of the scaling solu- 
tion, and is related to global properties of the gauge fix- 



ing, namely the Kugo-Ojima confinement criterion and 
the Gribov-Zwanziger horizon condition, see also [Io| . 

The functional DSEs and FRGs are derived from the 
effective action of the theory, expanded in its Green func- 
tions. With the abbreviation 
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the effective action is given as 
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with the gluon field A and the (anti-)ghost fields C, C, 
<fi = (A, C, C) and suppressed Lorentz- and colour- 
indices. In the third line in Eq.@ we have introduced 
an abbreviated notation which will be used throughout 
this work. The one-particle-irreducible Greens functions 
p(2n,m) w ith 2n external (anti-)ghost legs and m external 
gluon legs are the expansion coefficients of the effective 
action in the field expansion. It is convenient to repa- 
rameterise these expansion coefficients with the help of 
the coefficients 

p(2,0) = z (2,0) g(2,0) p(o,2) = £(0,2) g(j>,2) ^ 

of the kinetic terms. We then obtain the rescaled coeffi- 
cients r( 2n,m )(p) given by 



p(2n,m)^ _ p(2n,ro)^ 
2u 



(4) 



with P = (pi, ..., P 2n+m)- This parameterisation implies 
that the coefficients of the two point Green functions 

rim [p) = s ^\ P ), rM( P ) = s$*> (p)t (5) 

carry only the canonical momentum dependence of the 
kinetic terms and Z^°' 2 \ Z^ 2 ^ account for all quantum 
corrections. 

The reparametrisation (0J also entails that the Z- 
factors on the right hand side of @ carry the whole renor- 
malisation group scaling of the vertex functions r( 2n ' m ) 
in terms of the renormalisation scale fi. Together with 
the standard renormalisation group (RG) equation of the 
theory, 



r = o, 



(6) 
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we then learn from ([4]) and ([6]) that the expansion coef- 
ficients f ( 2 ™' m ) do not depend on fj,, i.e. 



djji 



0. 



Vn, m £ Z . 



(7) 



Note in this context that ^ is also valid in the pres- 
ence of the RG-adaptcd regulator terms in the functional 
renormalisation group equations (FRGs) 0, [l7| ■ This 
allows us to derive constraints for the vertex functions 

f (2n,m) also m the 

presence of the regulator, see [6j, lS| 
for details. In the present work we also show, see Ap- 
pendix [5] that the FRG-analysis of the IR-asymptotics 
can be further simplified, allowing for a more direct ap- 
proach. 

Eq. ([7]) already suggests that the RG-invariant coeffi- 
cients 

r (2n,m) do nQt 

carry any (global) anomalous scal- 
ing in terms of momenta. We shall show in the following 
that this is indeed the case. 



III. UNIQUENESS OF THE GLOBAL SCALING 



Then the scaling relations for the kinetic terms, the ghost 
and gluon dressing functions, read 



lim Z^°\\p) = \ K ^Z^°\p), 



\->0 



(11) 



The total global scaling A 2 ' 2 " " 1 of the full vertices 
p(2n,,ro) a j g0 mvo i V es the anomalous dimensions of the 
propagators and reads 

t2 n ,m = dm+m + ^ {^ nK 2fi + ni^Oa) + ^2n,m ■ (12) 

Previous analyses in 0, d, S EH were initiated similarly. 
However, instead of evaluating the deviation K2n.m to the 
standard anomalous scaling, the deviation K2n,m from the 
canonical scaling in the critical dimension four of Yang- 
Mills theory was evaluated, 



^2n,m — d>2n,va\d—A ~\~ ^2n,r 



(13) 



where 



For this proof we are only interested in the global scal- 
ing behaviour for the coefficient functions T^ 2n ' m \ Mod- 
ulo logarithms this entails the global scaling 



limr< 2 "' m )(Ap) = A 2(d2 



^£ n ' m) (p), (8) 



where r as stands for the infrared leading term, and p = 
(pi,...,P2n+m)- The coefficient 



d d-2 



(9) 



is the canonical scaling dimension of the vertex r^ 2 ™'" 1 ' 
with I = 2n+m. Note that it is only sensitive to the total 
number of external legs. It can be directly derived from 
Eq. (f3J): plugged into Eq. (T2J) it matches the canonical 
scaling of the (5-function of total momentum conservation, 
—d/2, and that of ghost, anti-ghost and gluon fields in 
position space, (d— 2)/4. Thus it matches the momentum 
scaling of the momentum integral f and that of 

the fields in momentum space. As a result the scaling 
§§§ includes the canonical momentum scaling of the one- 
particle irreducible Green functions as well as the scaling 
of the couplings. Hence, only in the critical dimension of 
Yang-Mills theory, d = 4, the canonical scaling dimension 
(|9"|) agrees with the classical momentum scaling. 

It turns out that the present parametcrisation (jSJ) in 
terms of K2n,m enables us to significantly simplify the 
proof given in Ref. @] of the uniqueness of the «2n,m- 
At its core the reason is the natural book-keeping of the 
necessary RG-scaling by the Z-factors that incorporate 
one factor of l/2«2,o or 1/2ko,2 for each external leg of 
the vertex r^ 2 ™'™). 



As for the most basic 



(cp. Eq. ©) 



we obtain by definition 



«2,0 = K .2 = 0. 



(10) 



Ad 2n +ro + 2 ( 2nK 2 : o + mK 0t2 ) + K 2n , m , (14) 



with the deviation Ad 2n +rn of the canonical scaling from 
that in the critical dimension four, 



Ad t =d l -d;| d=4 = (2-1)- 



(15) 



and I = 2n + m. This is adapted such that the K2n.m for 
primitively divergent vertices describe the full scaling of 
the corresponding dressing functions in d dimensions. 

We emphasise that in principle additional logarithmic 
scalings should be included into (JH) , (jTTJ) . However, even 
if present, additional logarithmic scalings do not change 
the relations between the 4,m and are therefore irrele- 
vant for the purpose of the present investigation. We also 
add that self-consistent logarithmic scaling laws have to 
satisfy additional consistency conditions. 

In four dimensions an explicit solution of the whole 
tower of DSEs and FRGs in terms of K 2n ,m was first 
derived in 0] and then generalised to d dimensions in 0] . 
It reads 



1 



; K o,2 



K2r, 



K 2,0 



(n — m)n + (1 — n) ( — — 2 



(16) 



In [8| we already gave a proof for the uniqueness of 
Eq. (|16[) . In the following we reformulate this proof in a, 
to our mind, more elegant and instructive manner that 
also allows us to reduce the number of presuppositions. 
We first derive constraints for K 2 „ jm from the functional 
renormalisation group and Dyson-Schwinger equations 
and then show that K2n,m — for all n, m. The resulting 
expression for the K2n,m from Eq. (|14p then agrees with 

Eq. mi). 
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FIG. 3: Functional renormalisation group equations for the gluon and ghost propagator. Filled circles denote dressed propa- 
gators and empty circles denote dressed vertex functions. Crosses indicate insertions of the infrared cutoff function. Only one 
possible insertion of the infrared cutoff function per diagram is shown. 



A. Constraints from the functional RG 

The FRGs for the ghost and gluon propagators are 
given diagrammatically in Fig. [3] With a mode cut-off, 
which only removes a single momentum mode, the regu- 
lator insertion is proportional to a 5-function and simply 
restricts the loop integral to a given momentum p 2 which 
we take to be vanishing. Then the loops on the rhs of 
the flow simply count the powers of global momentum 
scaling of the quantum fluctuations, no initial condition, 
similar to the classical term in the DSE, appears, see also 
Q. The potential cancellations necessary for the initial 
condition are discussed at the end of our proof. 

We are now counting anomalous dimensions on both 
sides of the equations in terms of powers of one external 
momentum scale p 2 in the infrared region p 2 <C ^qcd- 
For the global scaling ([5]) considered here all anomalous 
dimensions in terms of internal momenta of the loops 
translate directly into anomalous dimensions of the ex- 
ternal momentum scale. This is also true for the vertex 
equations considered below. In this respect the regulator 
insertion, denoted by the crosses, carry the anomalous 
dimensions of inverse propagators 0, Q . 

The constraint equations for K2 n ,m can be derived in 
several ways. A somewhat pedestrian approach is to 
count anomalous dimensions «2 n ,m of the dressing func- 
tions on both sides of the equations and then converting 
to K2n.m with the help of ([14"]) . More efficiently, we note 
that the K2,Oj k o,2 carry the renormalisation group scal- 
ing of the corresponding Green functions and match on 
both sides of the FRG equations. In particular this is 
true for the propagator FRGs in Fig. [3l Consequently 
all K2,0) ^0,2 drop out of the FRG-relations for a general 
vertex T2 n ,m- Note also that the sum of the canonical di- 
mensions d,2ij, in a given diagram for T2n.n1 simply 
gives the total canonical dimension (fen+m , and hence the 
d^ij also drop out of the FRG-relations. Then, we are 
left with relations for solely the K2%,j ■ For the propagators 
the constraints read 

= k ,2 = min (2« 2 ,i , 2k , 3 , k , 4 , «2,z) , (17) 
= K2,o = mm(2rt 2 ,i , «2,2 , «4,o) , (18) 



from the gluon and ghost FRGs. For the lhs of these 
equations we used that ^2,0 = ^0,2 = by definition, cp. 
Eg .pup. The minimum prescription on the right hand 
side of (|17p. (fTS]) takes into account that only one of the 
diagrams may be leading in the infrared. The constraint 
([15]) from the ghost-FRG entails 



«2,i > 0, 



K 2 ,2 > . 



S 4 ,o>0, (19) 



and at least one of these has to be zero for (fT5|) to be 
satisfied, 



«2,1 



0. 



or 



K 2 ,2 = . 



K 4 , = . (20) 



The same analysis 

K2,l, Ko,3, «0,4, ^2,2 > 



for (jTT)) entails that 
with at least one of them 
being zero. For the proof below, however, Eqs. (|19p and 
(|2U|) will be sufficient. 

We conclude the FRG-analysis with a discussion of the 
FRG relations for general Green functions. Schematically 
these relations read 



H2r. 



mffi (ft2„+2,m ' K 2n : m+2 , 



) 



(21) 



where the first two terms are the tadpole contribu- 
tion with ghost tadpole (K2n+2,m)i an d a gluon tadpole 
(«2n,m+2) respectively. The dots stand for other dia- 
grams with at least two vertices. It follows that K 2 n,m 
appears as the tadpole contribution in the relation for 
K2n-2,m and K 2 „, m _2, and more generally 



«2(r 



i-2s S K 2n ■ 



(22) 



for all r < n and 2s < to. This allows us to relate general 
K2m, m to either «2,i for odd m and 2s = to — 1, or «2,2 
for even to and 2s = m — 2. Thus we have 



> 



k-2,2 for to even 
R2 1 for to odd 



and we conclude with ([19]) that 

K2n,m > , Vn, TO S IN . 



(23) 



(24) 



in general space-time dimension d. This constraint to- 
gether with Eq. (|20p will be important in what follows, 
as it summarises in a closed form the infinite number of 
constraints from higher diagrams. 
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FIG. 4: Dyson-Schwinger equations for the gluon and ghost propagator. Filled circles denote dressed propagators and empty 
circles denote dressed vertex functions. 



B. Constraints from Dyson-Schwinger equations 

The Dyson-Schwinger equations for the ghost and 
gluon propagators are given diagrammatically in Fig. 31 
whereas the corresponding equations for the ghost-gluon 
vertex are displayed in Fig. O For the ghost gluon ver- 
tex we have two DSEs which are derived from either the 
functional gluon DSE or the functional ghost DSE, see 
Q. As already mentioned, the potential cancellations 
necessary for the classical terms are discussed at the end 
of our proof. 

We have seen in the analysis of the FRGs that the 
^2n.m-constraints boil down to simply counting the ver- 
tices involved in a given diagram and summing up the 
corresponding R2n,m- The same would apply to the DSEs 
if we only had dressed vertices in the DSE diagrams. 
However, there is always one bare vertex which then 
counts as K2n,m — K2n.m = —A^2n,m- These differences 
are given by 

AK 2 „, m = Ad 2 n+m + X ( 2nK 2,0 + TOK 0j2 ) . (25) 

For example, we are thus led to re^i — Are2,i for the right 
hand side of the ghost propagator DSE, and zero on the 
left hand side similar to the ghost FRG. This simple 
counting applies to all the diagrams. For its chief im- 
portance in the proof we introduce the abbreviation 



Ak = Are 2 i- 



(26) 



The constraints derived from the propagator DSEs dis- 
played in Fig. [4] are then given by 







K 2 1 — Aft, Kq,3 — Ako.; 



-Are , 



k oa - Ak oa , 2k 0j3 - Ak 0i4 J i (27) 
= mm(re 2 ,i - Are) . (28) 



Certainly, these relations can be derived as well in the 
pedestrian way of counting K2n,m on both sides of the 
equations and converting them to K2n. m - Note that in 
contradistinction to the FRG equations the DSEs do de- 
pend on K2.0 and ko, 2 via the Are2 n ,m- 

In the two different DSEs for the ghost gluon vertex 
in Fig. [5] we again apply the now familiar counting and 
obtain 

«2,i = mui (S 2 ,i + S , 3 - Ak , 2k 2 ,i - Ak , R 2 ,2 - Are) , 

(29) 

from the upper relation in Fig. [5] and 
K2.1 = min^2K 2 ,i - Ak ,3 , 2k 2 ,i — Ak , 

K 2j2 - Ak ,3 , ^4,o — Ak , two-loop^) ,(30) 

from the lower relation in Fig. [5j These constraints will 
be used in the next subsection. We emphasise that the 
above relations are valid in arbitrary dimensions as the 
FRG relations derived in the last Section fill Al In con- 
tradistinction to the FRG-relations the DSE relations de- 
pend on the dimension via the Are 2 n m- 



Proof 



We now proceed to show that 

K2n,m = 



(31) 



is the only scaling solution of infrared Yang-Mills theory 
in Landau gauge. To this end we note that scaling in the 
ghost-DSE (23) implies that 



Ak = K2 i > 



(32) 



where we have used that the ghost-FRG entails re 2 ,i > 0. 
Furthermore we obtain the two constraints 



^4,o > Ak , 



K 2 . 2 > Ak , 



(33) 
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FIG. 5: Dyson-Schwinger equations for the ghost gluon vertex. Empty circles denote dressed vertex functions. All internal 
propagators are dressed; the corresponding filled circles have been omitted for clarity of the figures. One-loop diagrams with 
the same scaling behaviour are only shown once. The ellipses denotes the other one- and two-loop diagrams which are not 
needed for our analysis. 



from the two DSEs for the ghost gluon vertex, where we 
use < K2,i < «2,2 — Ak from (|29j) . and < K2,i < 
K4 — Ak from (f5U|) . However, the ghost FRG led to 
the constraint l[20]). Together with (J32J) and ([33]) this 
immediately leads to 



Ak = 0, 



(34) 



and therefore also ^2,1 = due to 

It remains to be shown that this implies that all of the 
K<2n.m need to be zero. To this end we resort to the FRG- 
equations. The FRG-equations for all K2n.m contain at 
least one diagram that solely depends on iV 2njm ghost- 
gluon vertices. This implies 



K 2n,m < ^2n,ra K 2.\ — . 

Together with ([M]) this leads to 



&2r, 



0. 



(36) 



Written in terms of our original anomalous dimensions 
K-2n.rn, see (|14p . this implies that 

K2n, m = Ad 2n+m + i (2nK 2 ,0 + TOK ,2) , (37) 



with 



Arf 2 n+™ = (2 - 2n - m) 



d-A 



(38) 



sec (fl4"|) and ([T51 respectively. Eq. (j3"?| and (|38[) represent 
the unique scaling solution of infrared Yang-Mills theory. 

However, we would like to emphasise that the values 
of the anomalous scalings K2,o, K o,2 of the propagators 
cannot be fixed by scaling arguments. They have to be 
computed by solving the corresponding FRG and DSE 
equations, see e.g. U 0, [13, d, O ■ Instead, there is one 
final piece of information which can be extracted from 
the scaling analysis, namely constraints on the values for 
K2, 0,^0, 2^ From (|27|) we deduce Akq.4 < 0. Using 



we have (4 — d) 1 2 + 2ko,2 < 0, and with (|3"4")) we then 
conclude that M, [11 



^2,0 



1 14- 
2^ 2 " ~2~ 



(39) 



Eq. (|39p is in accordance with the Gribov-Zwanzigcr 
confinement scenario @, 0] which predicts ghost- 
enhancement. 

Finally we discuss the fate of the classical terms with 
0. The propagator scaling stemming 



^0,2 



class 
v 2,0 



from the quantum fluctuations is summarised qualita- 
tively as 



14-d 

Kon > 

2,0 _ 2 4 



K0,2 < 



A-d 



(40) 



(35) With k 0j2 < k 



class 
0,2 



0, the classical contribution does 



not change the gluon scaling. In turn, for k 2j o > 0, the 
inverse ghost propagator would be dominated by its clas- 
sical part, which has to be cancelled if scaling applies. 
This is the adjustment of the horizon condition, see e.g. 
0, [l(| ■ Note also that the bound (j40"|) in principle also al- 
lows for K2,o < for d < 4. In practice, however, one finds 
H2.0 > 0, i.e. K2,o(d = 3) = 0.40, and K2,o(d = 2) = 1/5 
for a classical ghost-gluon vertex [l8| . For these values to 
drop below zero, the full ghost-gluon vertex would have 
to deviate drastically from the classical vertex. This is al- 
most excluded by lattice computations [2l|, [52] and DSE 
self-consistency checks d, [23| . 

Numerical solutions for the ghost and gluon DSE 
as well as the corresponding FRGs in agreement with 
Eq J37J) and |35]) have been given in 0, El, [H, H, HI . 
In [26j a truncation has been employed which effectively 
converts the DSEs into equations with dressed vertices 
resembling the structure of FRGs. This procedure al- 
lows for interesting numerical solutions. However, the 
vertices used are neither consistent with the (unique) in- 
frared scaling laws Eq. (|3"T)) and (f3"8"|) . nor with the stan- 
dard RG-scaling. At its core this is due to the fact that 
the system of FRG equation cannot be solved by in the 
IR as the gluonic vertices in (2|| are too singular. 
These comments complete our proof. 



2n + m 2n + (m-l) 2n + (m - 1) 2n + (m - 1) 

# = ! 

2ra + m (2ro — 1) + m 

FIG. 6: General Dyson-Schwinger equations for p( 2n > m ). The straight legs denote either ghost or gluon lines with multiplicity 
given below the diagram, respectively. 



IV. KINEMATICS 

We have derived the unique scaling solution ([57)) with 
the assumption of conformal scaling. As already argued 
in [|| and in the introduction, possible conformal invari- 
ance entails that no further constraint on k 2 ,o an d Ko ; 2 
apart from (|40[) can be derived from pure scaling argu- 
ments. Here we shall make this formal argument more ex- 
plicit, and also discuss the question of existence of a spe- 
cific scaling solution including various kinematical limits. 

This is also interesting for the following reason: From 
our analysis in the previous section we know that the clas- 
sical purely gluonic vertices are sub-dominant, and the 
ghost-gluon vertex has its classical global scaling. Note 
however that a ghost gluon vertex with T^- 2 ^ ~ S^ 2 ' 1 ) 
requires the cancellation of the classical term in the DSE 
and the initial condition in the FRG equations in some 
kinematical limits, i.e. vanishing ghost or anti-ghost mo- 
mentum. In these limits the diagrammatical scaling 
would then be dominated by the classical part/initial 
condition. 

As already discussed above, the FRG equations can- 
not provide constraints on k,2.o, k o.2, as the k's drop out 
of the FRG-scaling analysis after employing the param- 
eterisation The FRG-diagrams then reduce to one 
loop diagrams with bare propagators and dressed ver- 
tices f ( 2rl ' m ) with only canonical global scaling due to 
([36]) . Here and in the following we use the term 'global 
scaling' when all external momenta scale in contrast to 
'kinematical scaling' when only some external momenta 
are involved 1 . Note that the fH 2 ™'™) carry kinematical 
singularities as is well-known from perturbation theory. 
In Landau gauge these may be meliorated by transver- 
sality. Since we are working in Landau gauge we are only 
interested in transversal contributions to all vertices; in 
Landau gauge the set of fully transversal vertices together 
with the ghost propagator and transversal gluon propa- 
gator provide a closed set of FRG and DSE equations. In 
turn, vertices with at least one longitudinal gluon satisfy 



DSEs and FRGs that also depend on the transversal ver- 
tices and only on the transversal gluon propagator, see 
also PH. 

For these reasons we concentrate on the purely 
transversal part of the vertices, T±^ 2n ' m '. For the treat- 
ment of kinematical scaling we also extend our notation 
to 

i Pr+l,2n i 2n+l,2n+s , P2n+s+l,2 n+m) 
A^O x 2(d r+a+K ;- m ) r^.™)^ t (41) 

where Fx^. 2 ™ stands for the infrared leading term, 
Pi j = [pi, ...,pj), and p is not exceptional by itself. Here, 
s counts scaling gluon momenta, and r — r± + T2 counts 
n scaling ghost momenta and T2 scaling anti-ghost mo- 
menta. Analogously to ([T^]) the total scaling A 2 ' 2 ™< m of 
the full vertex T 2 "' m reads 

*2n,m = d r+s + K 2n ,m ■ (42) 

The global scaling (|37[) is reproduced by 

K 2n,m = K 2n,m , V?' + S = 2n + 171 - 1 , (43) 

due to momentum conservation. 

In the following we shall show that a dressed ghost- 
gluon vertex without kinematic singularities is a possible 
solution of the DSE and FRG systems. In this case the 
scaling relations for the ghost-gluon vertex r^*- 2 ' 1 -' arc 

1,0 0,1 n { A A \ 

K 21 = k 2 \i = K 2,i = . (44) 
A ghost-gluon vertex with (|l4"l) reads 

r±£2c(p= ?) = K(p)^f(p + ?)r bc ■ (45) 

where 

n^(p) = v-^, (46) 



This corresponds to 'uniform' vs. 'collinear' singularities in the 
terminology of Il5ll . 



and / is a non-singular function of both momenta. Here 
p is the gluon momentum, q + p is the ghost momentum 
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(<Z+p)_ 



- - pi 



P2 



FIG. 7: Infrared leading diagram in the DSE for the ghost- 
gluon vertex. 



and q is the anti-ghost momentum respectively. Ghost- 
anti-ghost symmetry is implicit with 

n ^(p)^ = n ^(p)(p + q)m : f(p + q,q) = f(i,p + o) ■ 

(47) 

Now we assume that the leading infrared parts of all ver- 
tices are already provided by diagrams only depending 
on ghost-gluon vertices. This entails that, with the ex- 
ception of the ghost-gluon vertex, all ghost or anti-ghost 
legs of arbitrary vertices are proportional to linear pow- 
ers of the corresponding external momentum. The rea- 
son is that each of these legs is attached to an internal 
transversal gluon. With (|4"T|) we then can always rewrite 
this diagram as proportional to linear powers of the ex- 
ternal momenta. 

With this fact in mind we can deduce the scaling 
from the general form of the vertex-DSEs. Diagrammat- 
ically they are depicted in Fig. [6l We begin our anal- 
ysis with vertices where one external gluon momentum 
is vanishing. To this end we take the external momen- 
tum p = Pbarc °f the bare vertex to zero by keeping the 
other momenta at non-exceptional values. Extracting the 
linear momenta assigned to the internal ghost lines and 
projecting onto transverse components we get 



J dqq d - 1 J dfi,-^ 



I 



(1+k) (g+p)2(l+K) 
X 1v (<l+P)p<loZpai J ,i-» m -i(P,q,Pl ■ ■ 'Pan+m-l) >( 48 ) 

where J dfl q stands for the angular integration stemming 
from the ^-integration, and fj, = fj, m . 

As an example we discuss the simplest possible dia- 
gram for the case of no further external gluons, m = 1, 
shown in Fig. [7] with a similar momentum routing as 
in Eq. (|4"B"|) . In this diagram each of the three ghost- 
gluon vertices generates one four-momentum. Due to the 
transvcrsality of the involved two gluon lines these can be 
converted into the factor q v (q + p) p q a also appearing in 
Eq. (|48|) . For a non- vanishing integral we need a further 
power in q, otherwise the integral is proportional to p v 
and vanishes due to projection with 11^, (p). Indeed an 
extra factor pi ■ q for some i = 1, 2 is generated by the 
internal gluon. In the general expression (|48|) this extra 
factor is provided by the kernel 1. Counting powers of 
the scaling momenta p and q we find 



" l, 2n,l 



d 
2 



2k. 



(49) 



for m = 1. From (|48[) we have d/2 from the integration, 
—2 from the two denominators of the propagators, —2k 
from the two ghost dressing functions in the propagators 
and +2 from the four powers of q. The other momenta, 
Pi and P2 do not scale. The meaning of Eq. (|49|) is the fol- 
lowing: we only obtain an infrared dominated, divergent 
and therefore scaling integral for k > d/4, otherwise we 
cannot say something definite from the scaling analysis. 
Whether the integral is then vanishing or finite depends 
on the details of the kernel X and the angular integration; 
consequently the integral may or may not display kine- 
matical scaling for vanishing gluon momentum in this 
case. 

On the other hand, in the general case where further 
external gluons are present, i.e. m > 1, we have contrac- 
tions of q v {q + p) p q<T with I pap2 ... Pm (p,q,l 2 
These contractions generate terms proportional to p^ 
with i = 1, (m — 1) that do not vanish even when all 
external gluon legs are transversally projected. Thus we 
can have terms in the integrand proportional to q v p p q a - 
In order to decide whether the integral produces kincmat- 
ical singularities or not we count powers in scaling loop 
momenta, i.e. two factors of q, two ghost dressing func- 
tions, two denominators and the integration. We arrive 
at the condition 



d-2 



2k < 



K > 



d-2 



(50) 



Thus kincmatical singularities only occur for k > (d — 
2)/4. However, for the scaling of the vertex in terms 
of the external gluon momentum we have to count all 
powers of scaling momenta q and p and arrive at 2 



,0,1 



d-l 



2k. 



(51) 



for m > 1, in agreement with [1 51 ] for d = 4. 

From Eqs. (|49|) and (|5"Tj) we observe that the kincmati- 
cal gluonic singularities arc smaller or equal to the lowest 
possible global singularity -k as long &s k < (d — l)/2, 
i.e. 

d-l 

> -K -> k < (52) 



o.i 



Consequently the gluonic vertex dressing in the scatter- 
ing kernels in Fig.[S]does not lead to divergences for p — > 
in neither the gluonic nor the ghost DSEs. In particular 
this entails that 



1,0 



> 0. 



(53) 



2 Since we also count powers of p Pi this relation possibly involves 
the tensor structure of the vertex F2n,m- The only case where 
this is important and interesting is the special case of the three- 
gluon vertex To, 3. It is easy to see that the tensor structure of the 
bare three gluon vertex cannot have kincmatical singularities and 
therefore cannot take part in kinematical scaling. Since Eq. 1151 I I ■ 
however, involves an explicit scaling momentum with external 
Lorentz-index it has to represent the kinematical scaling of a 
different tensor structure of the threc-gluon vertex. 
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We conclude our analysis with a brief discussion of 
higher kinematical singularities, which is worked out in 
more detail in Appendix [51 There we derive the scaling 
relation 



' \ ( 2-d \ 

I mm I — ft, — - h 2ft I else 

(54) 

To summarise: for k < (d— l)/2 we have obtained kine- 
matical divergencies that are small enough such that they 
cannot invalidate the global scaling relation ([37]). Note 
however that the bound for the existence of a Fourier 
transform of the ghost propagator is ft < (d — 2)/2, see 
also [U . If k exceeds (d— 2) /2, the ghost propagator in po- 
sition space cannot be understood anymore as tempered 
distributions which is a necessary condition for correla- 
tion functions in a local quantum field theory. Hence the 
above bound k < (d — l)/2 is not relevant for this physi- 
cally interesting case. We conclude that the relation ([54)) 
for the kinematical scaling together with the global scal- 
ing Eq. (|3"T|) is a possible solution of the tower of DSEs 
and FRGs. 

We do not want to further this discussion, in particular 
as the actual numerical -and physically sensible- solution 
satisfy all of the above bounds, that is ft nU m < (d — 2)/2. 
We merely would like to mention that even in the case 
/« > (d — l)/2 one cannot conclude that the system is ill- 
defined. It only means that such an anomalous scaling 
cannot be captured by the DSE tower of equations, where 
the single diagrams do not entail the full RG-scaling, and 
hence also do not have manifest scale invariance in the 
scaling region. A consistent solution to such a system 
necessarily requires non-trivial cancellations between dif- 
ferent diagrams. In the ghost propagator DSE, where 
such cancellations cannot occur, they are not necessary; 
the ghost-DSE is well-defined even for ft > (c? — 1) /2 in the 
case of full kinematical scaling of the ghost-gluon vertex. 

From our kinematical analysis we find that in the case 
of k > (d— 2)/2 all vertices have kinematic singularities. 
As stated above the DSEs cannot be used anymore. The 
FRGs, however, can straightaway be solved by = 0. 

This is also the reason why the above solution, Eq. ([54"| . 
for kinematical scaling is not unique in contradistinction 
to the global scaling relation Eq. ([57)) . These comments 
complete our proof of existence. 

V. SUMMARY 

In this work we demonstrated that there is only one, 
unique global scaling solution of infrared Yang-Mills the- 
ory in Landau gauge. To this end we introduced a pa- 
rameterisation of the one-particle irreducible Green func- 



tions of the theory that splits off the renormalisation 
group (RG)-scaling and reduces the discussion to that 
of the scaling properties of RG-invariant quantities. To 
our mind this greatly simplifies the proof as compared to 
our previous work Ref. [8[. We wish to emphasise again, 
that the values of the anomalous scalings ft2,o 5 K o,2 of the 
propagators cannot be fixed by scaling arguments. These 
have to be calculated explicitly from the corresponding 
FRG and DSE equations. 

In addition we also demonstrated the existence of a 
specific scaling solution that includes 'kinematical scal- 
ing' in various kinematical limits. In contradistinction 
to the global scaling relation, kinematical scaling is not 
unique. 

In general, the method presented here is also ap- 
plicable to other theories. As an example we shortly 
discussed its applicability to scalar quantum field the- 
ories in the introduction; in 0] we discussed scaling in 
the gauge-Higgs theory. These applications may also be 
extended to Yang-Mills theories in other gauges. 
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APPENDIX A: INTEGRATED FLOW EQUATION 

The standard form in the flow equation is depicted in 
Fig [51 where t = Ink is the logarithmic infrared scale 




FIG. 8: Functional renormalisation group equation (FRG) 
for the effective action. Filled circles denote fully dressed 
field dependent propagators. Crosses denote the regulator 
insertion dtR. 

and the cross denotes the regulator insertion d t R- If Rk 
is chosen as a mode cut-off and simply removes one mo- 
mentum mode from the theory, the loop integrations in 
the FRG equations for vertices is reduced to a single loop 
momentum q 2 ~ k 2 and we can explicitly apply the in- 
frared power counting. This has been done in [8[. For the 
sake of the comparison with the DSE is also convenient 
to rewrite the FRG as follows [l(|, 
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Upon integration from an initial momentum scale k = A to k = this yields 

m = SM + QTrlnr( 2 >[0] + rcn) - ± £ ^Tr * - - d v T^ + rcn , (A2) 



where S A = (T A - Tr ln(r< 2 ' [0] + R) A - ren) entails the 
initial condition at k = A and the integral term on the 
rhs of (|A2|) is an RG-improvement term. If we perform a 
momentum rescaling as in © including that of A: A — > 
AA, the last term shows at most the same scaling as the 
first one. This entails up to renormalisation, that 



lim r\d>] 

X-»0 



-Trlnr (2) [(^] 



(A3) 



as far as infrared scaling is concerned. Taking a gluon 
or ghost derivative of (|A3[) leads to the diagrammatical 
representation of the infrared asymptotics of the flow in 
Fig. [J 



APPENDIX B: VERTICES WITH TWO OR 
MORE EXTERNAL SCALING LEGS 

In this appendix we investigate the kinematical scaling 
for vertices Y 2n ' m with s soft gluon lines, r\ soft ghost 
lines and r 2 soft anti-ghost lines with r = n + r 2 . Due 
to ghost-anti-ghost symmetry we can always choose r 2 > 
r\. We will explore some general situations and discuss 
exceptions at the end of this section. 

First we discuss diagrams with isolated external legs 
with vanishing momenta, i.e. diagrams where these legs 
are not neighbouring. These are simple. In the case of 
isolated external ghost lines no divergences are encoun- 
tered due to Eq. ([55)) . Diagrams with isolated external 
gluon lines scale like the corresponding diagrams with 
one vanishing momentum, i.e. Eq. (|49[) if no hard exter- 
nal gluon lines are present and Eq. (f5Tj) otherwise. This 
can be easily verified by going through some explicit ex- 
amples. In case some of the external gluonic momenta 
are parallel we always find Eq. (|5ip ; the argument is sim- 
ilar to that given below Eq. (f5Tj) . Having said this, we 
concentrate on vertices with only neighbouring legs with 
vanishing momenta for the remainder of this section. 

In the following argument we concentrate on the FRG 
diagrams which arc only one loop. A part of the related 
diagrams decay into three classes depicted in Fig. [5] The 
other diagrams involve vertices with soft and hard parts 
that cannot be separated in the above way. We will eval- 
uate these diagrams at the end of our discussion. The 
first two classes of diagrams in Fig. [5] summarise possi- 
ble diagrams of neighbouring vanishing legs for r 2 = t\ 
or r 2 — T\ + 2. The third class summarises those with 



T2 = Ti + 1. The total scaling of a vertex with n ghosts 

: I U 1 ; + : $3 &1 : ; 1 ; 



soft hard 



soft hard 



soft hard 



FIG. 9: Diagram classes with kinematical scaling of only part 
of the vertices. 



and r 2 anti-ghosts with r = ri + r 2 and s gluons is given 
by 

t r<s = d r+s + i {rn2.fi + sk . 2 ) , (Bl) 

where we have extended the scaling for ghost gluon ver- 
tices with r\ = T2 to non-existing vertices with r\ ^ r 2 . 
We are now in the position to deduce the scaling of ver- 
tices Y 2n,m with s + r\ + r 2 vanishing external momenta 
for gluons+ghosts+anti-ghosts respectively, by relating 
it to the global scaling (|F31[) of a vertex with only vanish- 
ing external momenta: we simply remove the hard parts 
of the diagrams and close the remaining open lines with 
some combinations of soft vertices and propagators, the 
scaling of which we then subtract. What is left is the 
scaling £ 2l J m of the diagrams in Fig. [9l Since this proce- 
dure is insensitive to the specific combination of vertices 
added, we simply take a minimal number of additional 
soft vertices for the explicit computations, without loss 
of generality. 

In the first diagram in Fig. [9] we substitute the hard 
part of the diagram by a soft full three gluon ver- 
tex, hence adding a further vanishing gluon momentum. 
From the resulting global scaling of £ r . s +i we have to 
subtract the global scaling to 3 of the three-gluon vertex 
leading to 



L 2n,m — <T,*+1 



0.3 



tr 



2k. 



(B2) 



In (|B2[) we have used that d r + s +i — = d r+s — 1. As al- 
ready pointed out above, e.g. substituting the hard part 
of the diagram with two ghost gluon vertices and one con- 
necting ghost propagator leads to the same result. The 
same is true for any other combination of vertices in the 
hard part of the diagram. The < in Eq. (|B2[) expresses 
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the fact that there are cases where the kinematical situ- 
ation of the same vertex may be represented by the first 
or second digram of Fig. [9] and either one may carry the 
leading kinematical singularity. 

To analyse the second diagram we have to distinguish 
two cases. For n = T2 we substitute the hard part by 
a soft full ghost gluon vertex, hence adding a further 
vanishing ghost momentum. Apart from subtracting its 
global scaling ti,\ — 1/2, we also have to take care of 
the kinematical intricacies of ghost lines. The hard part 
of the second class of diagrams in Fig. [5] also carries a 
momentum dependence that is linear in the momenta of 
the soft ghost legs. This has to be added separately, and 
we arrive at 

*2n,m — tr.s + 1 + 1 — *2,1 = t r>s — K , (B3) 

where we have used that d r + s +i = d r + s — (d — 2)/4. If 
ri = T2 — 2 the two connecting lines between the soft and 
the hard part of the diagram are two anti-ghosts and 
we have at least to substitute the hard part of the dia- 
gram with two ghost gluon vertices and one connecting 
gluon propagator. In our counting we again have to add 
two linear powers of momenta due to the two connecting 
ghost legs. The net result is the same as for the case 
ri = r%. Therefore the result Eq. (|B3[) with r = r\ + r-2 
summarises both cases. 

In the third diagram we substitute the hard part of the 
diagram by a soft full ghost gluon vertex, hence adding 
a further vanishing ghost momentum. Apart from sub- 
tracting its global scaling t2.i = 1/2, the hard part of the 
third class of diagrams in Fig. [5] also carries a momentum 
dependence that is linear in the momentum of the soft 
ghost leg which has to be added separately. We are led 
to 

tin,™ = tr+l,s + \ - *2,1 = t r>s - + . (B4) 

Note that t^n m = ^r+i,s with r + 1 = In. 

We summarise the results of this section with the scal- 
ing relation 

\-^T + ¥ r 2 = n + l, 

\ min (— k, ^ + 2k) else 

(B5) 



where we have used ([B2 ]) .([B3 ]) and (|B4)l . Note that one 
can swap between the different classes of diagrams by iso- 
lating one or two soft ghost or anti-ghost vertices within 
the hard part of the diagrams. This effectively removes 
these vertices from the counting. However, since the scal- 
ing relations is different for the three classes it could pay 
off in a more singular behaviour. It is easy to check that 
this is not the case, and hence (|B5[) represents the full 
maximal scaling of the diagram classes depicted in Fig. [9] 



We close with a discussion of the remaining diagram 
classes. A interesting specific case is depicted in Fig. [101 
This seems to entail that one can have neighbouring ghost 
vertices or anti-ghost vertices. However, one can show 
that this only is possible at the expense of additional 
loops as is present in Fig. [TUJ At its core this relates to 
the fact that the ghost number of vertices r 2,l,Tn vanishes. 
The scaling of the related sub-diagram is always positive, 
in the present case it is 1/2. Restricting ourselves again 
to the case r\ < T2 , we deduce that for T2 > r\ + 2 we are 
effectively reduced to the soft scaling of a diagram with 
r\ ghost and r\ +2 anti-ghost legs and hence Eq. (|B5|) ap- 
plies. Because of ghost-anti-ghost symmetry this covers 
the general case. 




soft hard 



FIG. 10: Example from an exceptional class of diagrams with 
kinematical scaling of only part of the vertices. 
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